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Abstract 

The geometrical structure known as the Tulczyjew triple has proved to be very useful in describ- 
ing mechanical systems, even those with singular Lagrangians or subject to constraints. Starting 
from basic concepts of variational calculus, we construct the Tulczyjew triple for first-order Field 
Theory. The important feature of our approach is that we do not postulate ad hoc the ingredi- 
ents of the theory, but obtain them as unavoidable consequences of the variational calculus. This 
picture of Field Theory is covariant and complete, containing not only the Lagrangian formalism 
and Euler-Lagrange equations but also the phase space, the phase dynamics and the Hamiltonian 
formalism. Since the configuration space turns out to be an affine bundle, we have to use affine 
geometry, in particular the notion of the affine duality. In our formulation, the two maps a and /3 
which constitute the Tulczyjew triple are morphisms of double structures of affine- vector bundles. 
We discuss also the Legendre transformation, i.e. the transition between the Lagrangian and the 
Hamiltonian formulation of the first-order field theory. 
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1 Introduction 

Variational calculus is a natural language for describing statics of mechanical systems. All mathemati- 
cal objects that are used in statics have direct physical interpretations. Moreover, similar mathematical 
tools are widely used also in other theories, like dynamics of particles or field theories. In classical 
mechanics variational calculus was used first for deriving equations of motion of mechanical system, 
i.e. the Euler-Lagrange equations. 

In numerous works by W. M. Tulczyjew, for example in the book [3U] and papers [371 [58 , 34, 351155]. 
one may find another philosophy of using variational calculus in mechanics and field theories. This 
philosophy, especially the one leading to the construction called the Tulczyjew triple, has been recently 
recognized by many theoretical physicists and mathematicians. The main advantage of the approach 
developed by Tulczyjew and his collaborators is its generality. For example, using the Tulczyjew triple 
for autonomous mechanics we can derive the phase equations for systems with singular Lagrangians 
and understand properly the Hamiltonian description of such systems. One can even discuss systems 
with more general generating objects than just a Lagrangian function, e.g. systems described by 
family of Lagrangians or a Lagrangian function defined on a submanifold. For the details we refer to 

m 

Another advantage of the Tulczyjew's approach is its flexibility. Being based on well-defined 
general principles, it can be easily adapted to different settings. No wonder that there were many 
attempts to generalize the Tulczyjew triple to more general contexts of mechanics on algebroids or 
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different field theories (see e.g [51 [381 EH] ) • In our earlier paper [H] we started with constructing a toy 
model of the triple for field theory in the simplest topological situation. 

The purpose of this work is to construct the Tulczyjew triple for first-order field theory in a very 
general setting, i.e. in the case where fields are sections of some differential fibration with no additional 
structure assumed. The origins of the geometric structures we study lie in the rigorous formulation 
of the variational principle including boundary terms. We pay much attention to recognize physically 
important objects, like the phase space, phase dynamics, the Legendre map, Hamiltonians, etc. These 
issues are usually not well elaborated in the literature, as the classical field theory models use to 
concentrate on the Euler-Lagrange equations. Of course, we recover also the commonly accepted 
Euler-Lagrange equations, this time without requiring any regularity of the Lagrangian. 

Classical field theory is usually associated to the concept of multisymplectic structure. The litera- 
ture on the subject is very rich, so we mention only a few main papers. The multisymplectic approach 
appeared first in [37^ and [23] . Then, it was developed by Gotay, Isennberg, Marsden and others in 
[TQl [TTJ [T21 [13] - The original idea of the multisymplectic structure has been thoroughly investigated and 
developed by many authors, see e.g. papers by Carihena, Crampin, Ibort, Cantrijn, De Leon [2, 3 ( 4 
and Echeverria-Enriquez, Munoz-Lecanda [6] for general analysis of the multisymplectic structure and 
its application to the classical field theories, and by Forger, Paufler, Romer [71 [3], or Vankershaver, 
Cantrijn, De Leon [44] for the discussion of more detailed problems associated to the structure. An 
interesting discussion of the problem can be found also in the paper [21] . The Tulczyjew triple in 
the context of multisymplectic field theories appeared recently in |28j . A similar diagram, however 
with differences on the Hamiltonian side, one can find also in [9]. Another approach to field theory, 
based on differential forms on fibre bundles, is present in works by Krupkova and collaborators, e.g. 

E3|2S1|2Z|). 

Our approach to the Lagrangian and Hamiltonian formalism developed in the paper is different. We 
do not use directly the multisymplectic formalism, building instead the triple out of natural morphisms 
of double structures of affine-vector bundles. Also, we do not use the framework based on Klein's 
ideology and do not concentrate on the Euler-Lagrange equations nor regular Lagrangians, since the 
phase dynamics is for us the principal object. Using the affine geometry as a tool and following guide- 
lines of variational calculus we arrive to spaces and maps on the Hamiltonian side of the triple. The 
variational problem we start with determines uniquely the phase space together with its canonical 
structure which is different from the one in [9]. Moreover, the canonical structure of the phase space 
is not a multisymplectic form, but a family of symplectic forms on fibres over the base manifold with 
values in the space of forms on the base. The two structures are related, but not identical. As far as 
we know, a similar research is being done independently, for instance, by L. Vitagliano |45j . and E. 
Guzman. Just before submitting this paper we spotted a preprint by Campos, Guzman, and Marrero 
[5] dealing with similar questions. 

The starting point of our studies is a locally trivial fibration £ : E — > M over a manifold M 
of dimension m, whose sections represent fields, and the corresponding bundle J X E of first jets of 
sections playing the role of kinematic configurations. A Lagrangian is a map L : S 1 E — ► f2 TO , where 
Vl k :— /\ k T*M is the bundle of k- forms on M. The phase space of the theory turns out to be the 
bundle V = y*E ® E C*(^ m_1 ) denoted simply V = M*E ® E n™" 1 . Here, V* E is the dual of the 
vertical subbundle V_E in TE and £*(Sl m_1 ) is the pull-back bundle of fi TO_1 along the projection 
(:E->M. 

The Lagrangian side of the Tulczyjew triple is constituted by a map 

a: J X P -> V*J X E (g)ji B fi m 

being a morphism of double structures of affine-vector bundles associated with fibrations over V and 
S X E. The vertical derivative dL of the Lagrangian is a section of the bundle V*J 1 i? f2 m — >■ ^E 
and the (implicit) phase dynamics is defined as a subset T> of $ X V being the inverse image by a. of the 
image of dL, i.e., V = oT x (dL{$E)). 

Similarly, the Hamiltonian side of the triple is built on another morphism of affine-vector bundles, 
fibred also over V and J 1 E, 
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where ^E is the 'affine dual' of ^E, i.e. the bundle of affine maps from }\E to ft™ e) , and P^E is the 
affine phase bundle of the affine line bundle 9 : ^E — > V , an affine analog of the cotangent bundle. A 
Hamiltonian is a section of the bundle 9, i.e. H :V — > \ Y E, so its affine differential dH can be viewed 
as a map dH : V — > PP E and defines a phase dynamics T> — /3" 1 (dH (V)) . It is obvious that, being 
presented in a coordinate-free form, the whole theory is covariant. 

We would like to stress that all the geometrical objects we construct are not just postulated ad hoc, 
but discovered by starting from natural general principles and rigorous investigations of the geometry 
which arises in this way. Another nice feature of our approach is that it admits a straightforward 
generalization for field theories of higher orders. This issues, however, we postpone to a separate 
paper. 

This work is organized as follows. In section [5] we present the conceptual background of variational 
theories which justifies mathematical constructions and physical interpretations in particular examples. 
After introducing some notation in section 02 we pass to the Lagrangian side of the triple in section |4l 
Section [5] is devoted to the Hamiltonian side of the triple. 

2 Variational calculus in Statics and Mechanics 

Let us start with the simplest case of statics of a mechanical system. We shall assume that the set 
of all possible configurations of the system is a differential manifold Q. The tangent and cotangent 
bundles 

T Q '■ TQ — > Q, n Q :T*Q^Q 

will also be used. In statics we are usually interested in equilibrium configurations of an isolated 
system, as well as a system with an interaction with other static systems. The system alone or in 
interaction is examined by preforming processes and calculating the cost of each process. We assume 
that all the processes are quasistatic, i.e. they are slow enough to produce no dynamical effects. Every 
process can be represented by a one-dimensional smooth oriented submanifold with boundary. It may 
happen that for some reasons, not all the processes are admissible, i.e. the system is constrained. All 
the information about the system is therefore encoded in three objects: the configuration manifold Q, 
the set of all admissible processes, and the cost function that assigns a real number to every process. 
The cost function should fulfill some additional conditions, e.g. it should be additive in the sense that 
if we break a process into two subprocesses, then the cost of the whole process should be equal to 
the sum of the costs of the two subprocesses. Usually we assume that the cost function is local, i.e. 
for each process it is an integral of a certain positively homogeneous function W on TQ. There are 
distinguished systems, called regular, for which all the processes are admissible and the function W is 
the differential of a certain function U : Q — > M. In this case U is called the internal energy function. 

An equilibrium point for the system is such a point q G Q that all the processes starting from q 
have positive cost, at least initially, i.e. for some sufficiently small subprocess with the same initial 
point. Usually we formulate only a first-order necessary criterion for the equilibrium point. It says 
that a point q is an equilibrium point of the system if 

W(Sq) > 

for all vectors Sq £ T q Q tangent to admissible processes. Vectors tangent to admissible processes 
are called admissible virtual displacements. The set of such vectors will be denoted with A. It may 
happen that AnTQ does not project on the whole Q. We have then the set of admissible configurations 
C = tq(A n TQ). For regular systems the equilibrium condition assumes the form 

dU(q) = 0. 

We examine the interaction between two systems by creating composed systems. We can compose 
systems that have the same configuration space Q. The composite system is described by the inter- 
section of the sets of admissible processes and the sum of the cost functions. We describe our system 
by making a list of all systems that, composed with our system, have certain admissible configuration 
q as an equilibrium point. We observe that at each admissible q all the external systems interacting 
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with our systems can be classified according to their influence on our system. Moreover, in every 
class we can find a regular system, therefore the whole class can be represented by the differential of 
the internal energy of that regular system. We call a force the class of external systems interacting 
with our system. The force is represented by a covector, i.e. an element of T*Q. Instead of making 
a list of all external systems in equilibrium with our system at the point g, we can give a subset of 
T*Q representing those systems. The subset of T*Q of all forces in equilibrium with our system at all 
admissible points we call the constitutive set. For a large class of static systems the constitutive set 
contains the complete information about the system. The passage from the triple (Q, A C TQ, W) 
describing our system to the constitutive set O C T*Q is called the Fenchel-Legendre transformation. 

Let us now use the above concepts to describe the autonomous dynamics of a non-relativistic parti- 
cle. For simplicity, we will consider only the unconstrained case. Let us assume that the set of positions 
of the particle (possible configurations) is a smooth manifold Q. There are at least two approaches to 
the problem. The first deals with the finite time interval [to, ti], while the second with the infinitesimal 
time interval represented by the Dirac (5-distribution at t. The finite case provides a useful represen- 
tation of objects coming from statics, while infinitesimal approach leads to differential equations for 
phase trajectories that are commonly used in physics. We skip the details of the construction and 
provide here only a summary of the results of both approaches. 

For the finite time interval the configuration space is the space of all motions, i.e. pieces of smooth 
curves parameterized by the time, 7 : [tojii] 3 t Q. This space is not a standard manifold any 
more, therefore we have to precise the notions of a smooth function, a tangent vector and a covector. 
We will work with functions (usually called functionals) of the form 

(2.1) S(j) = t L(t 7 (t))df. 

In the above formula, L is a smooth function on TQ, called the Lagrangian, and t 7 denotes the tangent 
prolongation of the curve 7. A curve in the configuration space always comes from a homotopy, i.e. 
from a smooth map 

X :l 2 3(M)^x(M)e(3- 

Restricting the domain of t to [to,ti] for every s, we obtain a curve in the space of motions that is 
smooth by definition. The choice is justified by the fact that the composition of the curve with any 
function of the form (|2.1I) is a real function smooth in the usual sense. 

A vector tangent to a manifold is usually defined as an equivalence class of curves. In our situation 
we can adopt the same definition. Working with equivalence classes is difficult, therefore we observe 
that each equivalence class at a configuration 7 can be conveniently represented as a curve 

(2.2) £7 : [t ,ti] — > TQ 

such that tq o (5 7 — 7. In differential geometry we define covcctors as equivalence classes of functions. 
Equivalence classes are again too abstract objects, therefore we need a convenient representation for 
covectors. The idea of such a representation is given by performing variation of the functional S and 
separating boundary terms like in the procedure of deriving the Euler-Lagrange equations: 

(2.3) (ss, 5 7 ) = r {£ £(t 2 7(t)), *y(*)>d* + (7>L(t7(*i), *y(*i)> - (PL(fir(to), *y(*o)>, 

J to 

where EL denotes the Euler-Lagrange variation of L that depends on the second prolongation t 2 7 of 
the motion 7, and VL is a vertical differential of L with respect to the projection tq. We see that the 
convenient representation of a covector is a triple (f,Po,Pi), where 

(2.4) / : [t , ti] -> T*Q, tt q o / = 7, p Q e T 7(to) Q, Pl G T 7(tl) Q. 
The evaluation between £7 and (f,po,pi) is given by 



(2-5) {{(f,Po,Pi),Sj)) = f \f(t),6j(t))dt + (pi,«7(*i)> ~ <Po,*7(<o)>. 

J to 
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The symbol ((•, •)) will denote the pairing between convenient representations. 

The elements /, po, and pi have physical interpretation. The curve / is an external force acting 
on the particle during its motion, po and p\ are the initial and the final momenta. The constitutive 
set consists of all triples (f,Po,Pi) such that the particle moves along the curve 7 = ttq o /, starting 
at 7 (io) with the initial momentum po and arriving to 7 (ii) with the final momentum pi, while it is a 
subject to the force / along the motion. Having the constitutive set, we can discuss isolated system, 
i.e. systems with the external force / = 0, as well as the system interacting with external forces of 
different kinds. We may interpret the momenta as a result of an interaction between the system and 
its past and its future. Therefore it makes no sense to keep the momenta equal to zero. The space of 
momenta is usually called the phase space of the system. We see that in our example the phase space 
is T*Q. 

The constitutive set, as defined above, is a complicated object. We would like to describe it in a 
more convenient way, e.g. using differential equations for curves in forces and momenta such that their 
solutions restricted to any interval [to, *i] he in the constitutive set. We can obtain such equations 
using the infinitesimal approach to the dynamics. 

Passing to the infinitesimal formulation, we replace the finite domain of the integration, [£o,ii]> 
with the Dirac's (5-distribution at the point t. We see that the configurations are now elements of 
TQ. Since the configuration space is again a manifold, we have natural notions of smooth functions, 
curves, tangent vectors, and covectors. The "internal energy" function is now just the Lagrangian, 
and its differential is an element of T*TQ. Virtual displacements are vectors tangent to curves in TQ, 
i.e. elements of TTQ. We observe, however, that studying convenient representations of vectors and 
covectors for the finite formulation gives interesting results also in the infinitesimal limit. A virtual 
displacement of a configuration t 7 (t) is a vector <St 7 (i) in TTQ such that 

(2.6) r TQ ((5t 7 (t)) = t 7 (t), TT Q (5t 7 (*)) = Sj(t). 

In turn, the infinitesimal limit of the virtual displacement Sj that we had for the finite time interval 
is an element tSj(t) <E TTQ such that 

(2.7) T TQ (tS 7 (t)) = Sj(t), Tr Q (t(5 7 (t)) = t 7 (t). 

The virtual displacement and its convenient representation are two elements of TTQ related by the 
canonical flip 

km ■ TTQ — > TTQ. 

The constitutive set for our regular system in infinitesimal setting is the graph of dL. Again, the 
convenient representation of elements of the cotangent bundle in the finite time interval formulation 
provides us with another useful interpretation of the constitutive set. For the finite time interval the 
evaluation of {f 7 Po,Pi) on Sj reads 

(2.8) f 1 (f(t), 5j)dt + (pi, 57(*i)> - <Po, *7(*o)>. 
In the infinitesimal case we get 

(2-9) </(*),*y(t)> + ^<p(*),*y(t)>. 

In the absence of external forces, another description of a constitutive set can be derived out of the 
equation 

(2.10) ±(p(t),S 7 (t)) = (dL(t 7 (t)),5t 7 (t)) 

The left-hand side is the so called tangent evaluation between a vector tangent to T*Q and a vector 
tangent to TQ such that they have common tangent projection on TQ. More precisely, if p : R — > T*Q 
and : R -» TQ are two curves covering the same curve 7 : R — > Q, then 



(2.11) 



((tp(t),t5 1 (t))) = ±(p(t),5 1 (t)). 
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Since the vectors tSj(t) and 5t"/(t) are related by the canonical flip kq, 

KQ (6t 1 (t)) = tS 1 (t), 

the differential of the Lagrangian dL(t 7 (i)) and the tangent vector tp(t) are related by the Tulczyjew 
olq (which is dual to Kq), 

a Q : TT*Q — > T*TQ. 

In this way we have obtained another description of the constitutive set, called the phase dynamics 
and given by the formula 

TJ*M D D = a Q 1 (dL{JQ)). 

If the system is autonomous, then the constitutive set for any time t is the same. The condition 
for a curve p : R D I — > T*Q to be a phase trajectory of the system is that 

Vi € I tp(t) e D. 

The dynamics D can be understood as a differential- algebraic equation for a pair of curves, / : R — > T*Q 
and p : R — >• T*Q, covering the same curve in Q. A curve 7 : / — >• Q satisfies, in turn, the corresponding 
Euler- Lagrange equation, if the curve I 3 t 1— > ag 1 (dL(7(i))) G TT*Q is the tangent prolongation of 
its projection to T*Q. 

External forces can be included in the picture as follows. Equation (j2.10j) completed with the force 
reads as 

(2.12) (f(t),Sj(t)) + l<p(t), a 7 (i)> = (dL(t 7 (t)), *t 7 (t)>. 

The set of all elements of TT*Q with fixed projections on T*Q and TQ is an affine space modeled on 
the appropriate fibre of T*Q. The force / can be therefore added to a vector tangent to the phase 
space. The map aQ can now be extended to the map 

5q : T*Q x Q TT*Q — > T*TQ, a Q (f,u) = a Q (u + /). 
The dynamics with external forces (see [33]) is a subset D of T*Q Xq TT*Q given by 

D = a Q - 1 {dL(TQ)). 

All the structures needed for generating the dynamics from a Lagrangian (without external forces) 
can be summarized in the following commutative diagram of vector bundle morphisms 

(2.13) TT*Q — ^T*TQ 




The map aQ is an isomorphism of double vector bundles. Recall that double vector bundles are, 
roughly speaking, manifolds equipped with two compatible vector bundle structures. The compatibility 
condition can be expressed e.g. as the commutation of the two Euler vector fields associated with 
these vector bundle structures. A precise definition of a double vector bundle together with its basic 
properties can be found in [301 1211 [19] . 
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The map oiq is also a symplectomorphism between the symplectic manifolds (TT*Q, Ajloq) and 
(T*TQ, wjq), where djLOQ is the complete lift of the canonical symplectic form ujq on T*Q and lojq 
is the canonical symplectic form on T*TQ. 

It may happen that the phase dynamics is an implicit differential equation, i.e. it is not the image 
of a vector field. In some cases, however, the phase dynamics is the image of a Hamiltonian vector 
field for some function H : T*Q — > M. So that we can write 



(2.14) 



V = Pq 1 {&H{T*Q)), 



where /3q is the canonical isomorphism between TT*Q and T*T*Q given by the canonical symplectic 
form ljq on T*Q, 

Pq '■ TT*Q — > T*T*Q, {P Q (v),w)=u Q (v,w). 
Let us recall that the canonical symplectic form ujq is defined as the differential 

(2.15) ujq = dd Q 
of the Liouville form -&q given by 

(2.16) #q{v) = {Tr.Q(v),T* Q (v)}. 

The structures needed for Hamiltonian mechanics can be presented in the following commutative 
diagram: 



(2.17) 



T*T*<3 



0Q 



TT*Q 



T*Q 





The map [3q is an isomorphism of double vector bundles. 

The formulation of the autonomous mechanics described above has at least two important features 
when compared with the ones in textbooks: it is very simple and can be easily generalized to more 
complicated cases including constraints, nonautonomous mechanics, and mechanics on algebroids [XT] 
QUE!]. And last but not least, we need no regularity conditions for the Lagrangian. The Lagrangian 
and Hamiltonian can be functions, but one of them (or both) can be replaced by families of functions 
generating Lagrangian submanifolds in T*TQ and T*T*Q, respectively. It happens e.g. in the case of 
a relativistic particle in the Minkowski space [42]. Moreover, the generating object for dynamics on 
the Lagrangian side can be replaced by a 1-form different from AL. It happens e.g. for systems with 
friction. The crucial role is played by two mappings: oiq and /3q. 

The two diagrams (|2.16|) and (|2.17|) glued together are called the Tulczyjew triple for mechanics. We 
would like to emphasize the fact that the Lagrangian side of the Tulczyjew triple is not postulated, but 
derived form the variational calculus. The Hamiltonian side, present only in infinitesimal formulation, 
comes from the fact that TT*Q is equipped with two Liouville structures, i.e. is isomorphic to two 
different cotangent bundles T*T*Q and T*TQ (see [43]). 

The paper is devoted to deriving the Tulczyjew triple for field theory, i.e. to the case where config- 
urations are sections of a certain fibration. Like in mechanics, the Lagrangian side of the triple appears 
as a result of the existence of so called convenient representations of equivalence classes representing 
tangent vectors and covectors. The Hamiltonian side is again related to the canonical isomorphism 
between certain cotangent bundles. 
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3 Notation 

The notation used in papers concerning the geometry of classical field theory is usually very compli- 
cated, because iterated tangent functors have to be used. In this section we will present a system of 
notation that will be used in the following sections. We will try to introduce some rules which make 
the notation systematic for the cost of the length of some symbols. 
For M being a smooth manifold of dimension m, we denote by 

r m '■ TM — > M and ttm '■ T*M — > M 

the tangent and cotangent bundles. Let U C M be a domain of local coordinate system (q l )™ =1 on M. 
We have the adapted systems of coordinates (q % ,qi) and (q l ,Pj) on t^{U) and (U), respectively. 
The manifold M is the manifold on which the field is defined, in applications it can be e.g. the space- 
time. We will assume in the following that the manifold M is oriented. This assumption allows us 
to use even forms instead of densities and odd forms (which are not commonly used). However, we 
will relax this assumption in the example to have clear physical interpretations of geometrical objects. 
Instead of the usual notation /\ fe T*M for the space of /c-covectors on M, we will use the shorter symbol 
Q k . The space Cl k is a vector bundle over M. The vector space over the point x 6 M will be denoted 
by VL X . For a local coordinate system (q l ) there is associated the volume form r\ = dq 1 Adq 2 A • • • Adq m . 
We will denote by r\i the contraction 1(^7)77. 

Let 

C : E -)■ M 

be a smooth locally trivial fibration with the total space of dimension m + n. The total space is the 
space of values of the field, i.e. a field is a local section of the bundle (. In applications the bundle 
£ can have additional structures. On an open subset V C E such that Q{V) — U we can introduce a 
local coordinate system (q l , y a ) adapted to the structure of the bundle. 

The space of vectors tangent to E and vertical with respect to the projection onto M will be 
denoted by ME. For the restriction of te to the space of vertical vectors we will use the symbol ve- 
The bundle 

v E :VE — > E 

is therefore a vector bundle. The adapted coordinate system on (V), coming from (q l ,y a ) on V, 
will be denoted by (q l ,y a ,Sy b ). We will need also the dual vector bundle 

p E ■■ y*E — > E 

with local coordinates (q l ,y a ,pb) defined on p^}(V). 

The space of first jets of sections of the bundle £ will be denoted J 1 £'. By definition, the first jet 
j 1 cr(e) of the section a at the point e is an equivalence class of sections having the same value e at the 
point x = ((e) and such that the spaces tangent to the graphs of the sections at the point e coincide. 
Therefore, there is a natural projection j 1 ^ from the space J 1 ^ onto the manifold E, 

j X C: ^E 3jV(e) t—t e G E. 

Moreover, every jet j 1 c(e) can be identified with the linear map from T X M to T e E being the tangent 
map Tcr restricted to the space T X M . It is easy to see that linear maps coming from jets at the point 
e form an affine subspace of the space Lin(T x M, T e E) of all linear maps from T X M to T e E. A map 
belongs to this subspace if, composed with T£, gives the identity. Using the tensorial representation, 
we can therefore write the inclusion 

i\E c T*M ® T e E. 

The affine subspace S\E is modeled on a vector subspace of maps having values in vectors tangent to 
E at the point e and vertical with respect to the fibration (. In tensor representation the model vector 
space is T*M ® V e E. Summarizing, the bundle S 1 E -> E is an affine subbundlc in the vector bundle 



C*(T*M) ® E TE — > E 
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modeled on the vector bundle 

C*(T*M) ® B V£ — > E. 

The symbol £*(T*M) denotes the pull-back of the cotangent bundle T*M along the projection £. In 
the following, we will omit the symbol of the pull-back, writing simply T*M (gig TE and T*M ® E V£\ 

Using the adapted coordinates (q t ,y a ) in V C E, we can construct the induced coordinate system 

"i) 



(q l ,y a ,y b j) on (J 1 C) 1 {V) such that, for any section a given by n functions a a (q l ), we have 



In the tensorial representation, the jet j 1 cr(e) can be written as 

dq ®W + W {q (x))dq ®W 

where we have used local bases of sections of T*M and TE coming from the chosen coordinates. 

In the following, we will have to use iterated tangent functors J 1 and V. All jet spaces we will use 
are spaces of jets of sections of bundles over M. All vertical tangent vectors are vertical with respect to 
the projection onto M. It means that ^VE is the space of jets of sections of the bundle (,ov E , therefore 
the projection onto VE will be denoted by j 1 (C o v E ). Similarly, V ' \ X E is the space of vectors tangent 
to S X E and vertical with respect to the projection onto M. The projection to S^E will be denoted 
by v$i E . Both spaces, VJ 1 !? and J 1 Vi5, play a very important role in the Lagrangian formulation of 
the classical field theory. We will discuss the structure of these spaces in the next section. Now, we 
would like to point out that the coordinate system {q 1 , y a ) on V <Z E allows us to construct coordinate 
systems: 

(q\y a ,y b „6y c ,6y d k ) on tft o v JlE )~\V) CV^E 

and 

(q i ,V a ,6y\y c j ,5y d k ) on {v E o j 1 ^ o v E )Y\V) C £\IE. 



4 Lagrangian formulation 

In the first-order field theory, a Lagrangian is a map from the space of first jets of sections of the 
bundle C to scalar densities on M, covering the identity on M. Since we assumed that M is oriented, 
we can identify densities with m-covectors. A Lagrangian L is therefore a map 

L : J 1 ^ — ► VL m 

covering the identity on M . The space S X E is often called the space of infinitesimal configurations for 
the first-order field theory. Let us recall that in statics and other variational theories all information 
about the system is contained in a constitutive set which is a subset of the cotangent bundle of the 
configuration space. In mechanics we developed another description of a constitutive set, using so called 
convenient representations of covectors. The convenient representation approach led to differential 
equations for phase trajectories of a system. In field theory we can adopt the same scheme. In 
the infinitesimal approach the space of infinitesimal configurations is a manifold and the role of the 
internal energy is played by the Lagrangian. Since variations are vectors vertical with respect to 
the projection onto M, the constitutive set is given as an image of the vertical differential of the 
Lagrangian. Taking into account that the Lagrangian has values in f2 m , we get that the constitutive 
set is a subset of \/*S 1 E ®^i E O m . In mechanics, studying convenient representations led us to the 
concept of momentum and external force. Let us do the same for the field theory. 

4.1 The phase space 

The first step is recognizing the phase space for the first-order field theory on the fibration (. We 
can use the calculus of variations as the guide-line for the problem. Let I? be a compact region in M 
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with the smooth boundary dD such that it is contained in a domain of coordinates U. We will do 
now the standard calculations in coordinates that leads to the Euler-Lagrange equations, but we will 
not assume that the variations vanish on the boundary. In the following, we will denote by (er a ) the 
functions defining local section er, and by (a a ,S<r a ) the functions defining its variation, i.e. a vertical 
vector field on E along the section a. The action functional S evaluated on a gives 



ID 

and the variation of S evaluated on the variation of a gives 

— So" M ^ 
v dy a dy b j dqi 

Using the Stokes theorem, we obtain 



(4-1) <JS,*7)= / (— + — — )<!■ 



(4.2, {SSM )=ixp°*-i-,-^y^L 



dyh 



The term integrated over D gives the definition of the external forces that are now interpreted as sources 
of a field, while the term integrated over the boundary dD gives the definition of the momenta. An 
object that can be integrated over the boundary of the region D is a (to — l)-form, therefore the results 
of the evaluation of momenta over the variations should lie in the bundle of (to — l)-covectors on the 
base manifold M. We point out, however, that the momenta should be evaluated on variations rather 
than on infinitesimal configurations. It is a special case of autonomous mechanics when variations and 
infinitesimal configurations are represented by the same geometrical object (tangent vector), therefore 
we can evaluate momenta on velocities. We cannot evaluate the momenta on first jets, at least not in 
the usual sense. The phase space for the first-order field theory on the bundle C is therefore the space 

(4.3) v = v*e (giEn 711 - 1 . 

It is a vector bundle over E. We will denote the corresponding fibration by 

7T : V -» E. 

Using a base of sections of the bundle pe and a base (rji) in f2 m_1 , we can construct local linear 
coordinates on 7r _1 (U), 

(4.4) (q\y a y b ), 
such that a section of the bundle n is represented as 

P j b(q,y)dy b ®rij. 

The role of external forces is played by the source of a field. The source is represented by a section 
of the bundle V*E (gi E n m -> M. 

All the above calculations can be done in a coordinate- free form. We postpone it to section l4~4l 

4.2 The structure of iterated bundles. 

Passing from (|4.1j) to (I4.2j) we have used implicitly a canonical mapping 

k : y^E — ► ^WE 



which is analogous to the canonical flip km : TTM — > TTM. The role of k in the field theory is similar 
to that of km in mechanics, where the map % is used to obtain a convenient representation of a 
variation of an infinitesimal configuration. The idea of using convenient representations comes from 
mechanics formulated for finite time intervals. In the field theory, instead of the finite time intervals, 
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we have the bounded domains of integration D C M. A configuration is then a section of £ restricted to 
D and its virtual displacement is represented by an equivalence class of curves in the space of sections. 
Curves in the space of sections come from vertical homotopies, i.e. maps 

such that D is contained in an open set U and J is a neighborhood of in ML The verticality means 
that for any t we have C(x( x > t)) — x - Fixing x we obtain a vertical curve in E x , while fixing t we obtain 
a local section of £. Restricting the domain of x 1° D C U, we obtain a curve in configurations. The 
curves are classified, as usual, with a use of some functions on configurations of the type of an action 
functional. We observe that equivalence classes are conveniently represented by vertical vector fields 
b<j along a sections a over D. In the infinitesimal approach, a configuration is the first jet ] l a(x) of a 
section and its variation is represented by a vector S'^a tangent to the space of first jets and vertical 
with respect to the projection on M. From the convenient representation we obtain the first jet of a 
vertical vector field along a section, i.e. an element of J 1 Vi?. 

Let us be more precise and define the map k using representatives of elements of S l \IE and VJ 1 ^. 
For a vertical homotopy 

X:UxI^E 

of local sections of (, and for a point xo £ M, we can create two objects. Taking the first jet of the 
section x i-> x(:c, i) & t x o, we obtain the curve 

1 1-> j x x(^o, *) 

in S X E which is vertical with respect to the projection onto M. The vector 

tj 1 x(a;o,0), 

tangent to this curve at t = 0, is an element of VJ 1 ^. On the other hand, we can first take vectors 
tangent to vertical curves t >-} at t = 0, obtaining a vertical vector field along the section 

x i ^ x{ x i 0)- The vector field 

x I— > tx(x-, 0) 

is a section of the bundle C ° Pe ■ VE 1 — > M. Taking the first jet of this section at the point xo, we 
obtain the element jHx^O; 0) of J 1 V_E. 

Definition 4.1. The map k : VJ 1 ^ — > ^ME is uniquely determined by 
(4.5) /s(tj l x(iBo,0)) =j 1 tx(a;o,0). 

The definition is correct as the both sides of (|4.5[) are independent on the choice of the representative 
X for an element of VJ 1 £'. 

Both spaces, J 1 V£' and VJ 1 ^, are double bundles in the sense that they carry the structure of 
two compatible fibrations. The space J 1 Vi? is fibrated over VE and the fibration is an affine bundle 
modeled on the vector bundle T*M ®\/ E VV_E — > VE 1 . The projection from J 1 \/E onto ME is j 1 (( o ve), 
since the first jets are calculated with respect to the projection on M . In the adapted coordinates we 
have 

(q i ,y a ,5y b ,y c j ,Sy d k ) ^ (q\y a ,5y b ). 

The space J 1 V£' is also fibred over ^E and the corresponding projection will be denoted J 1 ^. In the 
adapted coordinates the second projection reads as 

{q\y a ,8y\y c 3 ,8y d k )^{ q \y a iy c 3 ). 

The compatibility of the two bundle structures means that the model vector bundle T*M Cg>v_B VVE 1 is 
in fact a double vector bundle. For the concept of a double affine bundle and its model double vector 
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bundle we refer to |20) . We have therefore two commutative diagrams of bundle projections 



^VE 



(4.6) 



PE 




T*M ® VB WE 



VE 



T*M® E VE 




The space VS 1 E carries also the structure of a double bundle. One of the bundles is an affinc 
bundle and the other is a vector bundle. The vector bundle hbration 



V} X E — y S 1 E 



in the adapted coordinates reads as 



(q\y a ,y b J ,5y c ,Sy d k )^(q\y a ,y b J ). 
The second fibration is the affine bundle hbration 

VJ X C : VJ 1 ^ — y VE 

which we obtain applying the vertical tangent functor V to the projection J 1 ^ : S X E — y E. The model 
vector bundle for this affine bundle is 

V(VE ® E T*M) -y VE. 

In the adapted coordinates the affine projection VJ 1 ^ reads as 

(q i ,y a ,y b j ,5y c ,Sy d k ) (q\y a ,5y b ). 

The compatibility condition of the two projections can be expressed as the assumption that the model 
space V(VE ®e T*M) is a double vector bundle. The structure of V^E can be summarized in the 
following two commutative diagrams 

V(VE ® B T*M) 



(4.7) 



^E 




VE® E T*M 




where the second diagram represents the model double vector bundle. 

The map k defined earlier is a morphism of double bundles covering the identities on side bundles. 
On the level of the model double vector bundles it corresponds to the canonical flip n e ■ TTE —y TTE 
restricted to vertical vectors. 

Another example of a space equipped with the double structure of a vector-affine bundle is the 
space ] X V of first jets of the bundle £ o tt. As the bundle of jets it is fibred over V and the fibration is 
an affine fibration modeled on a vector fibration VV ®v T*M —> V. The vector bundle structure on 
\ X V comes from the jet prolongation of the vector bundle projection 7r : V — > E. The double bundle 
structure of $ X E can be summarized in the following two diagrams, 



(4.8) 




VP ® V J*M 




VE ® E T*M 
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where the second diagram represents the model double vector bundle. Elements of the model vector 
bundle VP ®v T*M can be added to elements of i^V. In VP ®v T*M there is a subbundle of vectors 
which are vertical with respect to the projection on E, i.e. vectors v such that V7r(w) = 0. Vectors 
vertical with respect to the projection on E have the first component tangent to the corresponding 
fibre of the bundle n. As usual, vectors tangent to a fibre of a vector bundle can be identified with 
elements of the fibre itself. Therefore, if \ln(v) — 0, then v can be identified with an element of 
V® E T*M = M*E®e ^ m_1 ®e T*M. Note that adding vectors vertical with respect to the projection 
on E does not change the right-hand side projection, i.e., if p is a local section of C ° tt, then 

J 1 jr(j 1 p(s)+«) = J 1 ?r(i 1 p(a)). 

In coordinates, if = (x l , y a ,p\, y c k,p l ds) and v = (at* , y a , v° bk ) i.e. v — v\ k dy b ® rjj ®dx k , then 

\ l p{x) + v = {x\y a ,p j b , y c k,p l ds + v l ds ). 



In section 14.11 we have observed that sources of fields are represented by sections of the bundle 
V*E® E tt m M. Since 

v*e ® E n m c y*E ® E vr- 1 ® E t*m 

we see that elements of the total space of the bundle of sources can be added to first jets from J : P 
without changing any of the projections. This operation will be needed in construction of the field 
phase equations with sources in section [ 



In the next section we will construct the main map of the Lagrangian formulation of the field theory 
that maps covectors on the space of infinitesimal configurations to their convenient representations. 
For that we will need an evaluation between the space J 1 P of first jets of sections of the bundle 

C o 7T : V — > M 

and the space S l \/E of first jets of vertical virtual displacements. More precisely, we will construct an 
evaluation between the bundle 

JV: ^T^^E 

and the bundle 

Pv E ■ ^VE — > S l E 
with values in the pull-back of the bundle il m — > M by C ° j X C- 

Let p : M D U — > P be a local section of the momentum bundle. We denote by a the underlying 
section of the bundle £, i.e. a : M D U — 5- E is such that p o tt = a. Let also 6a : M D U — > ME 
be a vertical vector field along the section a. There is a natural evaluation between V e _E and P e = 
V*£? <Ei S!" 1-1 with values in , therefore (p, 5a) is a (m — l)-form defined on U C M . W can define 
the evaluation between j 1 p(a;o) and j 1 Sa(xo) using the formula 

((j 1 p(^o),j 1 <5cr(a;o)» = d(p, 5a)(x ), 

so that the evaluation is a map 

<(.,.» :^Vxj 1E l l VE^n m . 
In coordinates, if a is given by local functions (er a ), if the momentum is represented as 

p\dy a ® 77, 

and 5 a as 



d 

5a"- 

then 



{p,5a)(x) =p\(x)5a a {x)r ll 
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and 



d(p.d(r) (.)',,) = ( ^j-(x )Sa b (x a ) +p j b (x )^j-(x ) } il- 



therefore the evaluation ((, )) in coordinates reads as 

(((q\ y a ,ph, V c kJ dm ), (q\ V a , Sy\ y%, Sy^))) = p l dl Sy d + phSy", 



4.3 The map a 

Let us now define the main geometrical object of the Lagrangian formulation of the first-order field 
theory. 

Definition 4.2. The relation 

a : ] X V — > y*i x E (g)ji E n m , 

given by the condition 

((u,k(w))) = (a(u),w) 

for all w having the same projection on S X E as u, will be called the Lagrangian relation. 
In this case the relation a is actually a mapping. In coordinates we have 

a : {q\y a ,p> b ,y c k ,p l dm ) .— > {q\ y a , y c k , ^p l dh p> b ). 

I 

The map a is a field-theoretical analog of the Tulczyjew cxm in mechanics. It relates covectors on 
the space of infinitesimal configurations which are elements of V*J 1 i? (8>ji£ il" 1 , to their convenient 
representations. In the simplest case, when sources of the field are equal to 0, a convenient representa- 
tion of a covector is the first jet of a section of the momentum bundle. If there are no constraints for 
infinitesimal configurations of the system described by the Lagrangian L, the constitutive set of the 
system is given as an image of J X E by the vertical differential dL (understood as a map from to 
V*J 1 i? (g)jE f2 m ). Using the map a we can obtain a convenient representation of the constitutive set 
as a differential inclusion which we understand as a condition for sections of the momentum bundle. 
This differential inclusion we will call the phase dynamics of the field. 

Definition 4.3. The phase dynamics of the field, when sources are equal to 0, is the subset V of S X V 
given by 

V = a-^dL^E)). 

The phase dynamics is also called the Lagrangian field equations. Let us note that obtaining the 
Lagrangian field equations from a Lagrangian is in our theory very simple. We do not require any 
regularity of the Lagrangian. 

Definition 4.4. We say that a section p : M — > V is a solution of the Lagrange field equations if 
In coordinates it means that 

Edp j b _ dL j dL 

~W ~ dy 1 ' P b ~ Wi 

The equations, known as the Euler-Lagrange equations for field theory, are consequences of the La- 
grange field equations. 
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The Lagrangian side of the Tulczyjew triple for the first-order classical field theory can be presented 
in the following diagram: 

(4.9) JI73 > VJ 1 ^ Q, m 




There is one projection in the above diagram that needs explanation. It is the projection 

Let us fix a point v € S 1 E and denote e = j 1 C( w ) an d a; = £(e). In the vector space \I V } X E there is a 
subspace of those vectors which are tangent to fibres of the projection j : £ : S 1 E — > E. Since each such 
fibre is an affine subspace of T*M <g> T e E, the space of vectors tangent to the fibre is isomorphic to 
its model vector space which is T*M <g) V e i?. An element of V* S 1 E ® fi™, treated as a linear function 
on My^E with values in fi™, can be restricted to the subspace of vectors tangent to the fibres. The 
restriction is an element of 

j x m ® y* e E ® r*™ » ® r*™- 1 = 7? e . 

Summarizing, the projection £ is a restriction of a covector to the subspace of vectors tangent to fibres 
of a certain projection. It provides the Legendre map, defined by the Lagrangian, from the space of 
infinitesimal configurations to the space of momenta, 

X-.^E^V, \(v)=Z(dL(v)). 

Sources of the field can also be included in the picture. Recall that sources are sections of the 
bundle V*E ®e ^ m M and that an element of the total space of this bundle can be added to the 
elements of J 1 "?. We define the extended map 5 : V*£ ® E fl m x v — f V*J 1 E ®ji E VT by the 
formula 

a(pjV) = "(p+jV)- 

Definition 4.5. The phase field dynamics with sources is the subset V of V*} X E ®^i E fl m Xp J 1 ? 
given by 

V = a- 1 (dL{i 1 E)). 

Definition 4.6. We say that a pair of sections p : M — > V and p : M ->■ \/*E <g> E Q m is a solution of 
phase field dynamics with sources if 

p{x) +} 1 p(x) e v. 

In coordinates it means that 

Edp>b _ dL ■ _ dL 
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4.4 Phase space: geometrical version 

In section I4.1[ using a coordinate calculation, we have split the differential of a Lagrangian into two 
parts: the Euler-Lagrange part and the total differential part. We have used the formula (|4.2|l to 
determine the phase space for the problem. Now we can do the same intrinsically, i.e. without using 
any specific choice of coordinates (see [4T|V 

Let us fix a point Xq in M. Any element v of VJ 1 ^ over Xq can be represented as v = k(j Sa(xo)) 
for some section da of the bundle C, ov e- The same can be done for any element of VJ'i? with the use 
of the isomorphism 

re (w) : S l VE — > VJ l E 

which is defined analogously to K, using representatives of elements of VJ'i? and j'Vi?. It is convenient 
to use the notation 5a l (xo) — 8<t(xo))^ Given a covector <p 6 T* Q M, we can choose a local 

function / on M such that f(xo) = and d/(xo) = f- Now we define 

(4.10) F(cp,v) = K (j 1 e(x )), where £(x) = f(x)Sa(x). 

It is clear that the value of F depends only on the covector <p and the vector v, and not on the 
representatives. We have then defined the map 

F : T* Q M x Vji^jJ 1 ^ — ► Vj.^jJ 1 ^ 

which is bilinear. In coordinates, if 

then 

F(<p,v) = <p j 8y°JL. 

We can see that the projection of i 7 '((/?, w) on \JE is zero, i.e. F(ip, v) is vertical with respect to j 1 ^. 

For a one-form /i on J 1 ^ with values in Q m , we define a one-form on J 1 ^ with values in J7 m_1 
by the formula 

(fi, F(<p,v)) = tp A (iFH,v). 

In coordinates, if 

/ x = (( M ) o d2/ a + (M 1 )adyf)®^ 

then 

The form is vertical, i.e. it vanishes on vectors vertical with respect to the projection j 1 ^. 

There is an operation of total differential dM defined for forms on jet bundles with values in fi fe . 
For example, if /i is a one-form on J l E with values in fl k , then djvfA 4 is a one-form on with values 

in ft k+1 given by the formula 

(d M fi(j l+1 a(x )),Sa l+1 (x )) =d((/iojV, fo'))(a;o). 

Applying djf to ifA*, for /i being a one-form on ^E with values in fi m , we get that d^FA* i s a 
one-form on J 2 i? with values in ft" 1 . Let 

= (r? - d A /iFM and P(p) = iir/z, 

where r 2 is the canonical projection 

if : J 2 £ — > ^E. 

We have 
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and both forms E(/i) and are vertical with respect to the projection on E. In particular, if 

v e MS 2 E and Trf (u) = 0, then 

(E(»),v) = 0. 

Indeed, let us take a representative Sa such that v = 8a 2 {xo). Since v is vertical, Sa(xn) = 0, so 

(4.11) (£( M ), v) = (/*,T7f («)) - (d M i F ,v) = 

(H,6a 1 (x )) - (d M i F , fc 2 (a;o)) = ( 5a 1 (xo) ) - d((ijrjti) ojV, fox 1 )^). 

The vector 5a 1 (xg) is also vertical with respect to the projection on E, therefore we can find a function 
/ vanishing at xq and a vector u — 5uj 1 (xq) such that 

5a\xo) = F(df(x ),u), 

i.e. we can write (for the first part of formula (|4.11[) ) 

(/j,,5a l (x )) = (fi,F(df(x ),Suj 1 (xo)) = df(x ) A (i F fi(j 1 a(x () )), tfw^zo)). 

Using the fact that /(xq) = 0, we can write that 

df{x ) A (i F fj.(j 1 a(x )), Su} 1 (x )) = d (/{ipfxo^a, Suj 1 )) (x ). 

Now, let us concentrate on the second part of formula (|4. 1 1[) . Since the form ip/j, is vertical, the value 
of 

(ipfiQ 1 a(x)),Sa 1 (x)) 

depends only on the jet ] 1 a{x) of the base section and on 6a(x). The value of the differential 

(4.12) d M ((^M°jV,5a 1 ))(x ) 

depends therefore on the second jet j 2 <r(xo) and the first jet Sa 1 (xq). This means that in formula (|4.12[) 
we can substitute the section Sa by the section x h4 £(x) — f(x)5uj(x) that covers the same section a 
and has the same first jet at xq. We can now continue the calculation started in (|4. 1 1|) : 

(4.13) ( M(j V(x )), So 1 fa) ) - d ((ipti o jV, fc 1 )) (x ) = 

d (f{i F fi ojV, fc 1 ) - (ipfio^a, kQ 1 ^) )) (x ) = d(i F( u ojV, /c^ 1 - kQ^Xxq). 

Let us now observe that f(x)5uj 1 (x) — k(J 1 ^(x)) is a vertical vector for any x, because f(x)Suj 1 (x) 
projects onto f(x)Suj(x) and kQ^x)) projects on 5oj(x). Using verticality of ipfi, we see that 

(iF/iOjV, fSu 1 - kQ 1 ^)) 

equals on the whole neighborhood of Xq, its differential is therefore equal to zero. 

We have shown that, for any one-form /x on S 1 E 1 the form E(n) is vertical. The form P(/x) is also 
vertical by definition. For \x — dL we can therefore define two maps: 

£(dL) : S 2 E — ► y*E ® E tt m , (S (L)(j 2 a(x)), Sa{x)) = ( E(dL)(fa(x)), Sa 2 (x)), 

and 

V(dL) : J 1 ^ — > V*E ® E n m ~\ {V{L){j 1 a{x)), *<r(x)) = (-B(dL)(j 1 a(a;)), ^(x)), 
such that formula (|4.2[) assumes the form 

(5S,5a)= f (£(dL)tfa(x)),Sa(x)) + f (T{dL)(] l a{x)), Sa{x)). 

J D JdD 
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4.5 Simple example: electrostatics 

Let us write the Lagrangian side of the Tulczyjew triple for Electrostatics. In this particular example 
we have included also sources of the field. Since it is a nonrelativistic theory, we have decided to use a 
three-dimensional affine space A as our playground. For the affine space tangent and cotangent bundles 
are trivial, so it is possible to present mathematical objects that appear in the theory in a simple way. 
The model vector space for the affine space A will be denoted with V. It is a three-dimensional vector 
space equipped with a symmetric non-degenerate positive definite bilinear two-form g representing the 
metric. For the purpose of this example we do not assume that A is oriented. For the affine space A 
we have trivial tangent and cotangent bundles, 

t a : A x V — > A, ir A : AxV* — > A. 

Moreover, because of the presence of the metric, there is a canonical isomorphism 

g:V^V*, „_►<,(„,.), 

and a canonical scalar density q £ V*. The symbol /\ V* denotes the space of odd k- forms on V. 

The potential of the electrostatics is a scalar field, therefore we take E = A x R and M = A. The 
main bundle of the theory is the trivial bundle 

pr A : A x R — > A. 

The space of the first jets of sections of the above bundle can be identified with 

]}E = Axlx V* . 

The first jet of a section ip at a point x is just (x, <p(x),d<p(x)). Since in our case fi m (M) — Ax /\^V* , 
we get that the Lagrangian is a map 

3 

L:AxE.xV* — > Ax f\ V* 
covering the identity on A which reads as 

L(x,r,(i) = 

If ip is a section of pr A (i.e. a function on A), then 

L(x,<p(x),d<p(x)) = ^{d<p(x),g- 1 (dip(x)))g. 
Let us look at the other spaces involved in the theory, namely 

3 3 3 2 

V* S X E ® O m (M) = AxRxV* x (R® f\ V*) x (V <8> f\ V*) ~ A x R x V* x f\ D V* x f\ Q V*. 
An element of the above space will be denoted by (a;, r, fx, T], i?). There are two projections: 

3 2 

v : A x R x V* x f\ V* x f\ V* — > A x R x V* 

and 

3 2 2 

^ : A x R x V* x f\ V* x f\ V* — > ixlx f\oV*. 
The graph of the differential of the Lagrangian is a subset defined as 

dL{A x R x V*) = {(x,r,fj,,0, .9 _1 (m)j0), x G A, r e R, n £ V*}. 
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We see from the above that the phase space of the theory is 

2 

P = AxRx f\ V*. 

We are looking for an equation for a section of the phase bundle over A, i.e. for a map 

(4.14) 49ii — > (x,(p(x),E(x)). 
The Legendre map that associates the phase to a configuration is 

2 

\ : A x R x V* — >4xlx /\oV*, X(x,r,p) = (x,r, g' 1 ^)^). 

Elements of f^V*, i- e °dd two-forms, can also be interpreted as vector densities. In the case of 
electrostatics, the two-form or the vector density that we obtain here, integrated over a surface in A, 
gives the flux of the electrostatic field through the surface. 
The space S X V is in our case 

2 2 

S X V = A x K x f\ V* x V* x (V* <8> f\ V*). 

An element of the above space will be denoted by (x, r,p, p,, v). Using the same symbols (x, r, p) does 
not lead to any confusion, because the objects denoted by those symbols are conserved by every map 
we use. The analog of the space of external forces for electrostatic field is 

3 

y*E <g> n m (M) = A x R x f\ V\ 
whose elements are (x, r, p). The extended map 5 reads as 

3 2 2 3 2 

A x K x f\ V* x f\ Q V* x V* x (V* <g) f\ V*) — > A x R x V* x /\ V* x /\ V*, 

(x,r,p,p,fj,,v) ' — > (x,r,n,Mt(u) - p,p). 
It means that the inverse image of dL(A x R x V*) by a is the set D of all (x, r, p,p, fi, v) such that 

(4.15) V = T\V)^ k\t(v)=p. 

The set D represents an equation for the section (|4. 14[) of the phase bundle: if r = <p{x) and p — E(x), 
we get 

H = d<p(x), v=] x E(x), 

and the equations are 

(4.16) E(x)=g-HMx))jQ, 

(4.17) dE(x) = p. 

It is easy to see that substituting ()4.16|) to (|4.17|) we obtain 

d(5 _1 (dv?(a;))jfl) = P 

which can be written as 

(4.18) (A^)fl = p, 

i.e the Poisson equation for the potential of electrostatic field produced by the charge density p. 

All the mathematical objects of the theory have clear physical meaning: the field itself is a potential 
for an electrostatic field. The electrostatic field is a covector field rather than a vector field, however in 
the presence of metric we can canonically translate one into the other. The phase is a vector density 
associated to the field and used to calculate the flux of the electrostatic field through a surface. 
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5 Hamiltonian formulation 

In this section we will construct the Hamiltonian side of the Tulczyjew triple. The name "Hamiltonian" 
is usually associated with the time evolution of the system. In our approach, the Hamiltonian side of 
the triple gives just another way of generating phase dynamics. The interpretation of the Hamiltonian 
itself strongly depends on the particular theory. 

The space of infinitesimal configurations S^E is an affine bundle over E, therefore it is necessary to 
use the affine geometry and the notion of affine duality. In the following we will recall this notion and 
construct an affine analog of the cotangent bundle. Since the affine bundle that appears in field theory, 
i.e. j 1 ^ : J X E — » E, has a reach internal structure, we decided to work first with simpler objects, and 
then apply the results to j 1 ^. 

5.1 The affine-dual bundle 

Let us first recall some facts from the geometry of affine spaces (for more details see e.g [H]). Let 
t : A — > N be an affine bundle modeled on a vector bundle v : V — > N. Let us also fix a one-dimensional 
vector space U . We will use the symbol v* for the projection V* ®nU — > N. 

The vector space of all affine maps from a fibre A q , q G N, to U will be denoted by Aff(A q , U). Every 
affine map has its linear part, therefore Aff{A q: U) is fibrated over Lin(V 9 , U) ~ V* ® U. Collecting the 
spaces Aff (A q , U) point by point in N, we obtain a vector bundle 

(5.1) i 4 : Aff (A, U) -4 TV 
and an affine bundle 

(5.2) 9 : Aff {A, U) — >V* ® N U. 

The fibration 9 is an affine bundle modeled on the trivial vector bundle 

pr 1 : V* <g> N U x U — >V* ® N U. 

In the case U = R, the space Aff (A q ,M.) is called the affine dual of A q and denoted A^. 

It is always useful to write geometrical objects in coordinates. We will use a set of coordinates 
adapted to the structure. Let (x l ) denotes a local system of coordinates in O C N. Choosing a local 
basis e = (e a ) of sections of V, the dual basis e = (e a ) of sections of V* , a reference section ao : O — > A, 
and a non-zero vector u G U , we can construct the adapted system of coordinates, 

(x\/ Q ) in A such that f a (a) = e a (a - ao (<?)), q = r(a), 

and the adapted system of coordinates 

(5.3) (x\ip a ,r) in Aff (A, U) 
such that 

ip (a) = ((f a e a (a - a (q)) + r)u, q = r(a) 

for ip G Aff (A, U). In coordinates, the projection 9 is expressed as the projection onto first two sets of 
coordinates (x z ,ip a ). 

In the family of all smooth maps from Aff (A, U) to U we distinguish maps $ : Aff (A, U) — > U 
which are affine along fibres of 9 and satisfy 

(5.4) &((p + u) = *(y>) - u. 

Property (I5.4j) implies that in every fibre 9~ 1 {p) of the fibration 9 there is exactly one point (p p such 
that ^(fp) — 0. It means that the set , 3/ _1 (0) is the graph of a section of the fibration 9. On the 
other hand, having a section E we can define 



¥ B fo>) = £(%)) - ip. 
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It is obvious that satisfies condition (|5.4p . Therefore we have a one-to-one correspondence between 
smooth maps which are affine along fibres and satisfy condition (|5.4I) on one hand, and smooth sections 
of the bundle 6 on the other. 

The differentials of maps satisfying property (15.41) are such covectors on Aff (A, U) with values in 
U that, restricted to vectors tangent to fibres of 9, they give —id. The submanifold of such covectors 
will be denoted K^id- It is a coisotropic submanifold of T*Aff (A, U) <8>Aff(A,l/) U with respect to the 
canonical symplectic structure on T*Aff(A, U) <E>Aff(A,c/) V with values in U. Using the local system of 
coordinates (|5.3[) , we can construct the adopted system of coordinates on T*Aff (A, U) ®mua,u) U : 

(x l ,ip a ,r, <Ti,f a ,p). 

In the above coordinates the submanifold K—id is given by the condition p = — 1 and the canonical 
symplectic structure on T*Aff(A, U) ®MUA,U) U reads as 

WAff( J 4,c/) = dcrj A da; 1 + d/ Q A dtp a + dp A dr. 

Having a coisotropic submanifold, we can perform a symplectic reduction. Leaves of characteristic 
foliation are orbits of the cotangent lift of the natural action of U on Aff (A, U). The reduced manifold 
will be denoted by PAff (A, U). In coordinates the reduction is the map 

K-u 3 (x\ip a ,r,a z J a ,-l) .— ► {x\<p a ,a u f a ) e PMf(A,U). 

Elements of PAff (A, U) can be interpreted also as equivalence classes of sections of the bundle with 
respect to the following equivalence relation. Since Aff (A, U) is fibrated over V* ®mU and the fibration 
is modeled on trivial fibration with the fibre being U, the difference £2 — ^i of two sections is a map 
from V* ®jv U to U. It is therefore clear what means that d(Y,2 — £i)(p) = for some p € V* ® U. 
We say that two pairs 

(pi, Si) and (p2, £2) are equivalent if and only if p\ = pi and d(£2 ~ £i)(pi) = 0. 

The equivalence class of {p, £) is sometimes denoted by d£(p) and called the differential of the section 
£ at the point p. The manifold PAff (A, U) is obviously fibrated over V* ®n U. The fibration is an 
affine bundle modeled on T*(V* ®n U) <E>Aff(A,!7) U — > V* ®n U. 

The above construction of PAff (A, U) is analogous to the construction of PZ for an affine bundle 
Z — > M, modeled on the trivial bundle M x R — > M given in [16]. The only difference is that M is 
replaced by the one-dimensional vector space U . 

The affine bundle PAff (A, U) is actually a double bundle. The second bundle structure is inherited 
from the double vector bundle T*Aff (A, U) <E>Aff(A,u) U. Let us first recall that the structure of the 
double vector bundle T*Aff(A, U) ®Aff(A,t/) U 

(5.5) T*Aff(A, U) ® mA ,u) U 

Mf(A,U) Mf(A,U)* ®jv U 




is given by the two commuting Euler vector fields associated with the two vector bundle structures: 
the canonical one, 

d d d 

and the second one, 
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The second projection <; can be understood as follows. The covector ip G T* Aff (A, U) ®U can be 
restricted to vectors tangent at ip to the fibre of the projection Aff (A, U) — > N, i.e. to the space 
Aff(^4 g , U). Any such fibre is a vector space, therefore vectors tangent to the fibre can be identified 
with elements of the fibre itself. This leads to the identification of the restriction of ip with an element 
of the dual to the fibre, i.e. an element of Aff(A q , U)* ® U. We denote by (r^)* the projection from 
Aff (A, U)* ® N U to N. 

The coisotropic submanifold K_id is an affine subbundle of the canonical bundle structure ^ Aff (A,u) 
and a vector subbundle (over a submanifold, see [19]) of the second bundle structure <j. The image 
<;(K_id) consists of all elements h of Aff(A q , U)* ®U that satisfy property (|5 .4|) : 

h(ip + u) = h{ip) — u. 

Elements of Aff(A g , U)* ®U satisfying property (|5.4p are in a one-to-one correspondence with elements 
of A q itself. A natural identification can be established as follows. Every a € A q gives rise to a linear 
map on Aff(^4 9 , U) with values in U by evaluation, i.e. 

h a : ip i — > tp(a), 

however h a does not satisfy (|5.4[) . Property (|5.4p is satisfied by — h a , i.e. the map ip h-> — (p(a). 

For dimensional reasons, every element of Aff(A qi U)* (g> U that satisfies property (|5 .4[) is of the 
form — h a for some a £ A q . 

The Euler vector field V2 is tangent to K^id and projectable with respect to the symplectic reduc- 
tion. It gives rise to a vector bundle structure 

PAff (A, U) -> A. 

We have therefore the double bundle 

(5.6) PAff(A,J7) 




N 



with the left projection being an affine bundle and the right projection being a vector bundle. 

We can repeat the above constructions for A = (j 1 C°C) _1 ( ;z; )i i- e - the fibre of the bundle S 1 E — > M 
that is fibred over = E x , V = ((op E )^ 1 (x), and U — il™. For simplicity, we will denote (j 1 C°C)~ 1 ( a; ) 
with {i^E). The space Aff(]\E, fi™) will be called the affine dual of i\E and denoted by i\E. Usually 
the affine dual of an affine space is the vector space off all affine functions on the affine space with real 
values. Here we replace real numbers with top forms on M, but we keep the name. As a result of the 
above general construction, we get the affine dual P X E together with the affine fibration 

6 X : J>. — > V, 

and the correspondence between sections of 9 X and maps which are affine along fibres of 9 X with values 
in f2™. Collecting affine dual spaces P X E point by point in M, we obtain a vector bundle 

j + C : J f £— 

and an affine bundle 

0: $E—>V. 



Sections of the fibration are in a one-to-one correspondence with maps from PE to f2 m covering 
the identity on M, affine along fibres of 8, and satisfying property (|5.4p . We have also the manifold 
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PJt E fibred over V and equipped with the canonical family of symplectic forms with values in Q m , 
parameterized by points in M and obtained by reduction from V* ^ E ®^ E f) m . 

We then expect that the Hamiltonian description of the first-order field theory will be connected 
with the fibration : S^E — > V and the space E. In particular, a Hamiltonian is a section of 
and the differential of this Hamiltonian at p £ V is an element of PJ^_E, therefore the image of the 
differential of a Hamiltonian is a subset of PJ^£\ Diagram (15.61) for S X E takes the form 

(5.7) pjt£ 




V S X E 




E 



Let us end this section by constructing coordinates in P^E adapted to the structure of the double 
bundle. In section [3] we have introduced coordinates (q* ,y a ,y j) in J 1 ^ and in section B~T1 coordinates 
(q l ,y a ,p 3 a ) in V . Since PE is fibred over V and the fibration is an affine bundle modeled on trivial 
bundle V x m ^ m — > V, it will be convenient to use coordinates (q l ,y a ,p J a , r) in ^ E such that r is an 
affine coordinate along the fibres of 9. If <p = (q l , y a 7 ph, r) and j 1 cr(x) = (q l , y a , y h j), then 

ip(j x a{x)) =p j by b 3 +r. 

In \/*PE Ojt_E ^™ we have therefore the adopted coordinate system (q l ,y a ,pi a ,r,£ a ,y b k, p) and the 
coisotropic submanifold K_id is given by the condition p = — 1. Diagram (|5.5[) in the case A = ^E 
takes the form 

(5.8) t/*Jt£:® Jt£; fi m 



that in coordinates reads as 
(5.9) 



.J f C 

E' 

(q\y a ,p> a ,r,Ca,y b k,p) 

"jtB / \ S 



(q\y a ,p i a,r) 



{q\y\y b 3 ,p) 



J+C (J f 0* 



(<f,y a ) 

After the reduction, we obtain coordinates (g l ,y a ,p J a , £ a ,y b k) in PJ^-E. Diagram (|5.7|) in coordinates 
takes the form 



(5.10) 



(q l ,y a ,P>a,£a,y b k) 

P8 / \ P? 



(q^y^p'a) 



{q\y a ,y b k) 



(q\y a ) 
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The family of symplectic forms with values in Q m parameterized by points of M is given in coordinates 
as 

(5.11) wpjt E = (d£a A dy a + dy\ A dp k b ) ® 77. 
5.2 The map /3 

In this section we will construct the map j3 which will be used in deriving the phase dynamics of the 
field from a Hamiltonian. It will be a field-theoretical version of /3q (see (|2.17[1 ). In mechanics there is 
a well-known formula relating Lagrangians to Hamiltonians, namely 

(5.12) H(p) = (p,v) - L(v). 

The origin of this formula lies in the procedure of composing symplectic relations (32j |33j [T] . Recall 
that a symplectic relation between symplectic manifolds (Pi, u>\), (P2, UI2) is a Lagrangian submanifold 
in (Pi x p2,wi — W2). If we deal with cotangent bundles, we can think of generating objects for 
symplectic relations. For example, it is well known that there is a canonical symplectomorphism 
TZf ■ T*F — > T*F* for any vector bundle r : F M . The graph of TZf is the Lagrangian submanifold 
in (T*P x T*F*,ujf — uip*) generated by the evaluation of covectors and vectors 

px m p* 3 (/,¥>)•— era 

in the following sense. The evaluation is a function defined on the submanifold F Xm F* C F x F* , 
therefore it generates a Lagrangian submanifold in T*(P x F*). The space T*(P x F*) can be naturally 
identified with (T*P x T*F*,ujf + wjr*). To get a Lagrangian submanifold with respect to the form 
uip we apply the transformation 

T*P x T*P* 3 (Ci, Ca) -> (Ci, -C2) e T*P x T*P 

to the generated submanifold. In the adapted coordinates (x l ,f a ) in F, (x 1 , f a , <n, ip a ) in T*P, and 
(x 1 , ip a ,(Ti,h a ) in T*P*, the isomorphism TZf reads 

■RF(x\f a , a u i/, a ) = (x\ t/, a , a,, ~f a ). 

If L is a function on F, the Lagrangian submanifold dL(F) C T*P can be treated as a symplectic 
relation between the cotangent bundle of the trivial one-point manifold and T*P. As a symplectic 
relation the submanifold is generated by —L. The operation of composition of symplectic relations 
does not in general lead to a symplectic relation, but since 72. is a diffeomorphism, we do not encounter 
such problems here. The set lZp{dL{F)) is a Lagrangian submanifold in T*P*. Composing symplectic 
relations means adding generating objects, we have therefore a generating family 

(5.13) H:F* x M F — > R, H(f,<p)=<p(f)-L(f), 

that in some cases can be reduced to a single function on F*. 

To clarify all sign problems let us mention that the family (|5.13[) is a generating object for a 
symplectic relation between a single point and T*P*. For a generating object of the Lagrangian 
submanifold 72._p(dL(P)) we have to take —H. 

In mechanics the whole procedure is applied to F = TQ. Then, — H is the generating object of 
a Lagrangian submanifold Dh in T*T*Q. The dynamics, i.e. a subset of TT*Q, is obtained as the 
inverse image of Dh by the map Pq : TT*Q — > T*T*Q, 

(5.14) Pq — o>q o IZjq . 
In this case /3q is also associated with ujq 



(5.15) 



Pq{v) = ujq{v, ■) for v E TT*Q. 
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In the classical field theory we need an afhne version of the above procedure, since in the space of 
infinitesimal configurations S^E we have only an affine structure on fibres over E. Moreover, we have 
to replace real- valued Lagrangians with Lagrangians taking values in the vector bundle f2 m . 

For simplicity, let us work first with an affine bundle A — >■ N and a vector space U like in the 
previous section. The canonical evaluation between elements of A and elements of Aff (A, U) is now a 
map defined on a submanifold A Xjy Aff (A, U) C A x Aff (A, U) with values in U. The evaluation in 
coordinates (chosen as in the previous section) reads 

Ax N Aff (A, U) 3 {x\f a ,<p a ,r) .— ► (f a <p a + r)u e U. 

The canonical evaluation generates a Lagrangian submanifold of 

T*(A x Aff (A, U)) ® ( AxAff(A,c/)) U 

with respect to the canonical U -valued symplectic form. There is an identification of the above cotan- 
gent bundle with 

T*A ® A Ux T*Aff(A, U) ® A ff(A,t/) U. 

The canonical symplectic form is identified with lua + w Aff (A, in- If we use the adapted coordinates 
(x l , f a 7 cTj, ip a ) in T*A®aU and appropriate coordinates (x n , (fi a ,r, cr-, h a , p) in T*Aff(A, ^)®Aff(A,i/)^i 
we get the generated Lagrangian submanifold given by the conditions: 

x l =x'\ a i = -a' i , ip a = (f a , f a = h a , p = l. 

To get a symplectic relation out of that submanifold, we have to change signs in the fibre of 
T*Aff(A, U) ®mua,u) U. The graph of Ha is therefore given by the conditions 

x = x , a-j = a,-, y a = tp a , J =—h, p = —l. 

We observe that the relation TZa is a map defined on the submanifold K_id with values in T*A (&a U. 
In coordinates, 

n A {x\ (p a ,r,ai,f a ,-l) = (x\-f a ,(j l ,ip a ). 
Since TZa is constant on leaves of the characteristic foliation of K-id, it reduces to a symplectomorphism 

K A ■ T*A ® A U — > PAff(A, U). 

For any function L : A — >■ U, from dL(A) C T*A <S>a U we obtain a Lagrangian submanifold of 
PAff(A, U). As a generating object we can choose a family of functions on Aff (A, U) parameterized by 
elements of A: 

H :Ax N Aff {A, U) — > U, H{a, ip) = L(a) - <p(a). 
Note that H satisfies condition (|5.4j) . i.e 

H{a, ip + u) — H(a, ip) — u. 

It generates, of course, a submanifold of T*Aff(A, U) ®U which, after reduction, equals TZA(dL(Aj). 
If the family reduces to a single function, the latter corresponds to a certain section of the bundle 
9 : Aff (A, U)-+ V* <g> U. 

Applying again the above constructions to A = J;j._E and recalling that we use the notation 
Aff(J^£ l , fi" 1 ) = JJ.-E, we get a diffeomorphism 

ftji^ : V*.) 1 ^ <x> fl m — > P^E 

which, restricted to every fibre over M, is a symplectomorphism with respect to appropriate symplectic 
f2™ valued forms. We will denote with (3 the composition 



(5.16) 



P = aoU^ A , /3:^V — > P^E. 
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In local coordinates, we get 

(5.17) (3(q\y a ,ph,y c k,P l ds) = (q\y a ,ph,Y,P h ck,y d i)- 

k 

The map /3 constitutes the Hamiltonian side of the Tulczyjew triple for the classical field theory. 

(5.18) Pjt£^ ] X V 




E > E 



The above construction shows that, in the first-order field theory, a Hamiltonian is not a density 
valued function on the phase space, but a section of certain affine bundle over phase space with one- 
dimensional fibres. Differentials of such sections are elements of an affine analog of the cotangent 
bundle. From the construction we obtain a family of Hamiltonian sections parameterized by elements 
ofJ 1 ^, 

Y> H : ^E x E T^ tfE, 
that corresponds to a family of density valued maps 

H : ^E x E J+£ Q m , Htfa, <p) = £(jV) - (ptfa). 

In some cases the above family reduces to a single generating section. It happens e.g. in Electrostatics 
(see section loT4|) . In such cases we obtain V from the image of the differential of the Hamiltonian 
section H : V — > S^E by means of the map f3. More precisely, 

(5.19) V = /T^dtffP)). 

Like in the Lagrangian case, the process of generating phase dynamics from a Hamiltonian is very 
simple. 

5.3 Structure of the phase space 

The phase space V is fibred over E. The fibration is a vector bundle. The space of vertical vectors \IV 
is therefore a double vector bundle fibrated over ME and V, 



VP 




V \JE 




E 



We define a one-form -Q-p on V with values in fi m 1 by the formula 

$ v (8p) = {v v {5p),Ve{8 v )). 
In coordinates, for Sp — (x 1 , y a , %P a , 8y c , 8p k d) , we get 

dv(8p) = p l a Sy a T] t , i.e. i? P = p\dy a (g> m- 
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The form fi-p is an analog of the canonical Liouville form on a cotangent bundle. Applying cim to "d-p, 
we obtain a one- form on J 1 ?' with values in Q m which in coordinates reads 

&W&-P = p\idy a ®r/ + p j b dy b j <8> rj 

and which is an analog of dj^M, the Liouville form on TT*M. 

Applying vertical differential d to (Im^p, we obtain a two-form on $ X V which can be treated as a 
family of presymplectic forms with values in fi m , parameterized by points of M: 

wji-p = dd.M'&v = dp l ai A dy a ® rj + dp J b A dy b j ® 77. 

It is easy to see in coordinates that 

f3*ujpjt E = —wji-p. 

There is an alternative construction of the map j3 that uses the language of differential forms on fiber 
bundles. The crucial role in the construction is played by the two-form di)-p with values in il m_1 [46 . 

5.4 Electrostatics 

The Hamiltonian side of the Tulczyjew triple for Electrostatics is simplified, because fields are sections 
of the trivial bundle ( : A x R — >• R. The bundle J 1 ^ is therefore a vector bundle: 

J 1 (:J 1 £ = 4xlxr — >4xl. 

The affine dual to this vector bundle is again the Cartesian product. Note that £l m ~ A x /\ 3 V* . We 
have the following identifications: 

S 1 E ~ixlx V*, 

2 3 
Jt£ ~ A x R x f\ V* x /\ V, 

2 

P~AxRx f\ V*. 
The bundle 6> : i ] E — > V is trivial, i.e. 

2 3 2 

9 : A x R x f\ V* x /\ F* — >ixlx A° y *' 

There is a natural action of /\ 3 T^* in the fibres of the above bundle by addition: 

(A 3 o^*) x [A x R x A 2 oV* x A 3 V*) — >ixlx/\ 2 V*x/\ 3 
(u, 0, r,p, A)) 1 — > {x, r,p,X + u). 

The canonical evaluation between (x, r, fi) G J 1 ^ and (n, r,f>, A) G J^-B with values in A 3 oV* reads 

(5.20) ((x,r,(i),{x,r,p,X)) =/iAp + A. 

The above evaluation generates a relation 1Z between the cotangent bundle of the space of infinitesimal 
configurations, 

3 2 

V*.! 1 ^®^" 1 ~ AxRxV* x f\ V* x f\ V* 
and the cotangent bundle of the affine dual, 

23 3 
V*^E ® Q m ~ 4xIx/\/'x/\ fxrx/\/'xl. 



28 



Katarzyna Grabowska 



An element (x, r, p, c, oS) <E V*J 1 i? ® 17 m is in the relation 7?. with an element (x, r,p, A, a, b, c) G 
V*jt£:®« m , if and only if 

p = cu, a = —a, b = p, c—1. 

There is an action of /\ 3 ^* m the cotangent bundle V*J t _E ® f2 m lifted from the action in J+.E. The 
image of 7?. (defined by c = 1) is invariant with respect to the lifted action and the quotient space 

2 3 

P^E can be identified with A x R x f\ V* x f\ Q V* x V*. 

The graph of 1Z is also invariant with respect to the lifted action, therefore there exists the quotient 
relation between V*J 1 i? <g) f2 m and P^E. This quotient relation is actually a map, 

TZ : y*^E^n m — > PJt£, 

(x, r, ,u, cr, uj) i — > (x, r, u, -a, p). 
Composing the map 7Z with a from the Lagrangian side, we get the map /3 

/3 ■. y*E (g) n m x E J 1 ^ — ► PJ f £, 

i.e. 

A x M x A 3 V"* x A 2 V"* x y x 7* ® A 2 oV" — > A x M x A 2 V* x A 3 V^* x V*, 
(x,r,p,p,p,u) i — ► (x,r,p,p- Alt(v),p). 

Equation (|4.17p for a section of the phase bundle V — > A can be generated from a Hamiltonian. The 
generating family 

h(x, r, p, X, p) = L(x,r, p) — p A p — \g 

reduces to a section H of the bundle 6>. Since the bundle is trivial, this section can be represented by 
the map 

3 1 
H:V^/\ V*, H{x,r,p) = -(*p)Ap. 

where * is the Hodge-star associated with the metric g. The Hamiltonian H generates the subset 

dH(V) C P^E, &H{T) = {(x,r,p,a,b) : a = Q,b = *p}. 

The inverse image of dH(V) by f3 is 

(5.21) p- 1 (dH(V)) = {(x,r,p,p,p,v); p = AU{v), p = *p}. 

Comparing (|5.2ip with (I4.15|) . we see that 

D = p- l {&H(V)). 



6 Conclusions 

We have constructed the Tulczyjew triple for the first-order field theory, starting from fundamental 
concepts of calculus of variations. Our results can be summarized in the following diagram of affinc 
and vector bundle morphisms 



(6.1) 



PJt£- 



Po 
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The left-hand side of the triple is Lagrangian, the right-hand side is Hamiltonian, and the phase 
dynamics lives in the middle. The phase dynamics V being a subset of S X V is interpreted as a condition 
for first jets of sections of the momentum bundle. It can be obtained from a Lagrangian as 

V = a- 1 {&L{S 1 E)) 

or from a Hamiltonian as 

P = /?- 1 (dE H (7')). 

The Hamiltonian is a section of an affine bundle 9 with fibres modeled on the vector spaces of volume 
elements on M. 

The Lagrangian side of the triple looks similarly in almost all papers devoted to the Tulczyjew 
triple for field theory. However, we would like to emphasize that all the spaces and maps that we use 
are constructed, not postulated, and have clear interpretation in the language of variational calculus. 
It is interesting also that the phase space V is not dual to the space of infinitesimal configurations ^E. 

To construct the Hamiltonian side of the triple we have used the notion of affine duality. Geometrical 
language for the first-order theory is not the only place in classical mathematical physics where affine 
structures are needed. Constructions that were used in section [5] are similar to those needed in time- 
dependent mechanics |16] or in the intrinsic formulation of Newtonian mechanics |15j . We expect also 
that the notion of affine duality will play an important role in higher-order theories. The Hamiltonian 
formulation of the first-order field theory is based on the canonical isomorphism PpE ~ V* S 1 E <S> fi m 
generated by the evaluation between S^E and its affine-dual bundle. 

In the Tulczyjew triple for mechanics, all three spaces, T*T*Q, TT*Q, and T*TQ are isomorphic. 
It is not the case in the field theory. The middle space is not isomorphic to P^E and V*J 1 E (gi f2 m . 
To have an isomorphism between all three spaces, we could replace the space S l V with the quotient 
space with respect to a certain equivalence relation. We decided not to do it, because passing to the 
quotient we would loose the obvious interpretation of the constitutive set as a first-order differential 
equation. However, it is clear that only very special differential inclusions T> C $ X V are generated by 
some Lagrangian or Hamiltonian. 

The spaces P$E and \J*S l E (g> VL m are equipped with canonical two forms that, restricted to every 
fibre, are symplectic forms with values in Q m . We have also a canonical presymplectic forms on fibres 
of the bundle J 1- P — > M. The phase space V possesses a canonical one-form with values in fi™ 1-1 which 
is an analog of the canonical Liouville form on T*M. 
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